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We review the current state of experimental knowledge on double cc production in e+e annihilation. The 



large cross-sections (O(20 fb)) for e" 1 



7* — > (cc) res processes have been confirmed by detailed tests 



and reproduced by a second group: they should now be considered well-established. The latest experimental 
results concern the case where the second cc system is above open-charm threshold: hidden-charm states continue 
to play a prominent role in the mass spectrum. Some "loose ends" in the field are also briefly discussed. 



1. Introduction 

Electron-positron annihilation to charmonium and 
an additional hidden-charm state, or a pair of open- 
charmed mesons — double cc production — has been 
known for only five years, but is now an established 
object of theoretical and experimental study. The field 
is still being driven by data. This presentation surveys 
current experimental knowledge: after reviewing the 
history (Section |2| and the now-established results on 
double-charmonium production e + e _ — > ip^ (cc) res 
(Section [3j, we discuss the latest experimental work, 
which focusses on states above open-charm threshold 
(Section [4|. Some "loose ends", which might reward 
renewed attention, are also noted (Section [5]). In clos- 
ing (Section |6| we summarize both established and 
new results, and their relation to theory. 
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Figure 1: The mass spectrum of the unreconstructed (re- 
coiling) system X in e + e~ — > J/tpX production in early 
Belle data |4j . The curves show a fit to a series of charmo- 
nium states (solid) over a continuous component (dashed); 
the shaded histogram shows the distribution for events in 
the J/tp mass sidebands. The charmonium lineshapes are 
broadened by initial state radiation: see the text. 



2. History 

This field grew from studies of inclusive charmo- 
nium production e+e" -> V W X [H|2]. An old CLEO 
analysis on a small sample [3] had presented evidence 
for direct decays T(45) — » J/tpX, distinguished by 
J/tp momenta above the endpoint for T(4S) — > BB[^ 
J/tp X] . Using data from early B-factory running — 
29.4ftT 1 on the T(45) resonance, and 3.0fb _1 in 
the continuum 60 MeV below — Belle p~] excluded such 
production, setting a limit B(T(4S) — » J/ipX) < 
1.9 x 10 -4 at 95% confidence. More importantly, they 
established picobarn cross-sections for e + e~ — ► ijy^ X 
processes in the continuum, and a peculiar momentum 
spectrum for the produced ipv>: in the J/tp sample, 
the cross-section fell to zero well below the momentum 
endpoint (see Fig. 3 of Ref. [1]). 

The c + c~ initial state being known, a simple rescal- 
ing of the J/tp momentum gives the mass of the system 
recoiling against it in the final state: 

M recoil (V) = yV*-^) 2 -^) 2 (1) 

An upper bound on momentum thus corresponds to a 
lower bound on the mass of the recoiling system. Stud- 
ied in this way [I], the surprising conclusion (Fig. [TJ 



was that the interaction e + e~ — > J/tp X does not pro- 
ceed if the mass of the recoiling system X is below cc 
threshold, while immediately above threshold there is 
significant two-body production, with X — rj c , Xco, Vc- 
Although familiar now, it is worth remembering how 
completely unexpected this result was at the time. 

It was then accepted [5] that charmonium produc- 
tion at y/s « 10.6 GeV was dominated by the e + e~ — * 
ipgg process, with a momentum spectrum extending 
to the kinematic endpoint and thus a M Teco n spectrum 
to low masses; e + e~ — > ipcc was an 0(10%) correc- 
tion. An additional contribution from e + e _ — > ip <?, 
where the charmonium develops from a colour-octet 
cc pair, was expected at momentum endpoint (low 
Mrecoil)- The contradiction with data led to reexami- 
nation of both theoretical and experimental methods. 

In the recoil mass technique, the system X is not re- 
constructed and identification is thus indirect. Alter- 
native explanations of the data have sprung from this 
limitation, while improvements to the method have fo- 
cussed on partial reconstruction and constraints. The 
original analysis [4! already used a mass-vertex con- 
straint in J/if> reconstruction, improving M rcco ;i res- 
olution by a factor of two. Contributions from QED 
processes are more troublesome. Initial-state radia- 
tion (ISR) leads to a high-Af roco ;i tail (typically model- 
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dependent) on the lineshape of any peak; more than 
four tracks are required to suppress contributions from 
other (low-multiplicity) QED interactions. 

The leading alternative interpretation of the data 
relied on such a process: e + e~ — > 7*7* t/il |]. 
While e + e~ ^7* — > tp X production requires the sec- 
ond state to be even under charge conjugation, £j£ = 
+ 1, the two- virtual-photon process allows ^ = ±1 
and in particular permits e + e~ — > 7*7* — > J/V 1 J/tp. 
As only the J/ tp rj c signal was statistically significant 
in the 2002 analysis [4 J , it was attractive to ascribe the 
peak (in part) to events with a second J ftp- A more 
ambitious (and essentially mirror-image) proposal [7] 
was that e + e _ — > tp gg did in fact dominate charmo- 
nium production, with the gluons sometimes coupling 
to a glueball state close to the rj c mass. 

These and other interpretations were effectively 
ruled out by an updated Belle analysis [5], using 
155 fb -1 of data and techniques designed to dis- 
criminate between the theoretical options. Three 
e + e _ — > tp rj c events were fully reconstructed, to be 
compared with the 2.6 ± 0.8 expected from the inclu- 
sive yield. The recoil mass scale was also calibrated 
using c + c~ — > Jisr tp' events and found to have a 
bias of less than 3MeV, ruling out confusion between 
rjc and J/tp. Attempts to include e + e~ — ► iptp^ 
components in the recoil mass fit resulted in negative 
yields, and restrictive upper limits on their contribu- 
tion (see Fig.[2j upper plot). An eventual confirmation 
by BaBar 9 found comparable results (lower plot). 

Belle also performed an angular analysis of J/tprj c '^ 
and J/ipXcO events. The distinctive forward peak in 
J/tp production angle (cos6* proc i — + 1) predicted by 
the two-virtual-photon model for J/tp J/tp (Fig. 2 of 
Ref. [5]) was not seen in data: results for all three 
states were consistent with shapes 1 + a cos 2 0, and 
equal coefficients for production- and helicity-anglc 
distributions (see Table [I]) as expected for a single vir- 
tual photon. The J/tp rjc fits were consistent with 
a = +1 (P-wave production) as required by rj c quan- 
tum numbers, and strongly disfavoured the —0.87 ex- 
pectation for a spin-zero glueball [7] . One- rather than 
two- virtual photon production (or the glueball expla- 
nation) is thus favoured by all experimental tests. 



Table I Coefficients from fits of the function 1 + a cos 6 
to J/tp production (# pr od) and helicity angle (#hei) data at 
Belle H], for e + e~ — > J/tp (cc) rea . See the text for results 
under the constraint ahci = a prod, and the expectation. 
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Figure 2: Baseline e + e~ — > J/tp (cc) r es results from Belle 
(0! upper plot), and the confirmation from BaBar ([9], 
lower plot). In the upper plot the dashed line shows the 
upper limit contribution from final states where the fitted 
yield is insignificant or negative: J/tp (J/tp, \ci,c2, tp'). 



3. Baseline results 

Belle [H] and BaBar [5] measurements of e + e~ — » 
ip('~* (cc) res cross-sections are summarised in Table [nj 
with competing explanations excluded, and reason- 
able agreement between the two experiments, these 
results are no longer in serious dispute. Remarkably, 
there seems to be no suppression of radially-excited 
states: cross-sections for tprj c , tp7j' c , tp' rj c , and tp' rj' c 
are all comparable. This presumably contains some 
hint as to the production mechanism. 

Low-order perturbative calculations, as embodied 
in nonrelativistic quantum mechanics (NRQCD, |101 
[TTJ[T2]) underestimate the cross-sections by an order 
of magnitude or more; the discrepancy is reduced, but 
not removed, when relativistic corrections are taken 
into account [TU]. Neither are other features of the 
data well-explained: NRQCD predicts a ~ +0.25 for 
e + e~ — > J hp XcO [TO] , disfavoured by the Belle analy- 
sis (Table p|, which prefers pure S-wave production. 

A calculation in the light-cone formalism [13 , how- 
ever, appears to match at least the J/tprj c cross- 
section. A variety of theoretical approaches are cur- 
rently being pursued, and it is no longer easy to char- 
acterise the issues at stake: on comparison of NRQCD 
and light-cone estimates, sec for example the very ex- 
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Table II Double-charmonium production cross-sections from Belle and BaBar, with theoretical pre- and post-dictions. 
Due to background-suppression criteria, the experiments report effective cross-sections for the case where the unrecon- 
structed state (ccjres decays to at least 2 ("> 0") or 4 ("> 2") charged tracks, and thus underestimate the cross-section. 
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tended discussion in Ref. [2] . Analysis of the merits 
of this theoretical work is outside the scope (and com- 
petence) of this review. 1 It does however seem clear 
that experimental work is still driving, rather than 
being directed by, theoretical studies. 



4. The new cutting edge: 

states above open-charm threshold 

Active experimental work has now shifted to the 
case where the system recoiling against the ip is above 
open-charm threshold. Even allowing for the various 
e+e" -> DWBW continuum components in fitting the 
inclusive M roco ii(^) spectrum, Belle found a 5.0a peak 
at (3936 ± 14) MeV (Fig. 1 of Ref. [IB]). The back- 
ground under this peak (called X(3940)) being too 
large for further detailed study, Belle explicitly recon- 
structed a D-meson and then considered cases where 
the remaining system was close to the D or D* in 
mass: M rcco ii('0 D) ~ m D (.>. Constraining such cases 
to match the mass also improved the M reco ;i(^) 
resolution. This allowed the reconstruction of a clear 
X(3940) — > D*D peak, and an upper limit on the same 
structure in DD (Fig. 3 of [16]) 

Belle has released an updated analysis [17] based on 
systematic use of this D(*)-tagging technique. Fig. [3] 
shows the M Teco n{ipD^) spectrum for 693 fb -1 of 
data, after reconstruction and mass-constraint of the 
J/tp, and then a D°, D+, or D*+ meson. A simulta- 
neous fit with the D(*)-mass sidebands is performed: 
clear and significant peaks are seen, corresponding to 
processes e + e~ — > ?/?DD, ?/>D*D, and -0D*D*. 



1 We note that a theoretical analysis has appeared since the 
workshop 15 , claiming that the discrepancy in the J/i/jrj c 
cross-section between NRQCD and experiment is now resolved. 
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Figure 3: The mass of the system recoiling against fully- 
reconstructed J/ip and D (upper plot) or D* + (lower plot) 
mesons at Belle [T7]. The shaded histogram shows the 
distribution in D^*'-mass sidebands: the fitted curve is 
discussed in the text. 



Monte Carlo study shows that we indeed expect 
these processes to be reconstructed in this way, with 
recoil mass resolution of about 30 MeV: smaller than 
the difference in D and D* masses. Disjoint samples 
|M r ecoii(^DW) -mt ag | < 70 MeV where the unre- 
constructed system is tagged as a D or a D* are thus 
selected (ISR leads to a 10% t/>DD ->■ -0DD* cross- 
feed), and M rGCO ii(-0 D^*-*) is then constrained to the 
mass of the tagged meson. This improves the resolu- 
tion on Af(DWl)M) by a factor of 3-10: results for 
the three samples are shown in Fig. [4] Peaks above the 
background are seen near threshold in each sample. 

Combinatorial backgrounds are taken into account 
via simultaneous fits to the data in the reconstructed 
D(*)-mass signal and sideband regions. The ex- 
cess over background is fitted with the sum of a 
threshold function to represent non-resonant e + e~ — * 
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Table III Resonant enhancements in e + e~ — > tjj D'^D'*' 
production at Belle [17]. The V an d the meson shown are 
fully reconstructed and then refitted under a mass con- 
straint: the remaining meson is identified as discussed in 
the text. 
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Figure 4: M(D ( * ) D ( * ) ) spectra for e+e" -» J/^D^dW 
events at Belle [17], where the first D-meson is explicitly 
reconstructed, and the associated D-meson is recovered 
via recoil mass selections and constraints: (upper) DD, 
(middle) DD*, and (lower) D*D* samples are shown. Dis- 
tributions in the reconstructed D'^-mass sidebands are 
shown in yellow, and reflections (middle) DD — > DD* and 
(lower) DD* — * D*D* in green. The solid curves show 
the fits described in the text, and the dashed curves their 
combinatorial background and reflection components. The 
second panel in the upper and middle plots shows the data 
and fit with these components subtracted. 



^,jj(*)E)(*) production, and an S-wave relativistic 
Breit-Wigner. In all cases the threshold term is in- 
significant, and a significant Breit-Wigner peak is 
seen. In A/(DD) the peak is broad, and the fit un- 
stable against variations of its conditions; in M(DD*) 
and A/(D*D*) the fit is stable and the resonant peak 
significant at over 5a. Various cross-checks are per- 
formed, including Monte Carlo and data tests of the 
background shape in D(*)-mass sideband and signal 
regions; study of charged- and neutral-D subsamples; 
and fitting of events with reconstructed D* , and asso- 
ciated D: the latter obtains results consistent with the 
J/?/>DD* analysis of Fig. [I] (middle), but with lower 
efficiency and significance. Parameters for the reso- 
nant enhancements are summarised in Table [TTTT 

The X(4160) enhancement has not previously been 
reported. The X(3940) mass and yield results are 
consistent with those of the earlier analysis [16] , while 
the width is larger than the published value of (15.1 ± 
10.1) MeV: a likelihood function non-parabolic in the 
width parameter had been noted in that case, with 
a 52 MeV upper limit at 90% confidence. The corre- 
sponding limit in the new analysis is T < 76 MeV. 

Cross-sections for J/^X(3940) and J/V>X(4160) 
production (Table III last row) are in the 20 femto- 
barn class, as for all the significant e + e~ — > ijj^> (cc) res 
processes seen to date. 



5. Sidelines 

The concentration on experimentally fruitful prob- 
lems establishing and measuring quasi-two-body 



processes e T e — » ipW ( cc )r 
neglect of inclusive e + e~ — 



:s — has led to a relative 
ip^X studies. The 2002 



Belle analysis [4 established the fraction 
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Figure 5: Inclusive e + e~ — > J/t/j X data, showing the com- 
binatorial background (yellow) and non-J/ip (cc) re3 con- 
tributions (blue) below open-charm threshold. (Adapted 
from a figure by P. Pakhlov.) 



Those results establish that e + e~ — > ipcc domi- 
nates charmonium production in the continuum at 
sfs sa 10.58 GeV, while tpw (cc) res cross-sections are 
at the 20 fb level, with no suppression of radially ex- 
cited states. Recent work shows that prominent reso- 
nant contributions continue above open-charm thresh- 
old, with similar cross-sections; this process is proving 
fruitful in the search for new hidden-charm states. 

The ball now lies in the court of theory. Interpre- 
tive work on the new states Y(3940) and Y(4160) is 
already underway, but a predictive account of e + e^ — > 
X amplitudes is still lacking. An advantage of the 
NRQCD approach is its pretension to universal appli- 
cation: at the Tevatron, for example, and the LHC. 
We await an accurate account of double-charmonium 
production in e + e _ annihilation that can also embrace 
quarkonium production at other facilities. 



using the following method: reconstruction and mass- 
constraint of J/ip, an d then an associated D^*- 1 - 
meson; rejection of contamination from c + c~ — > BB 
events using momentum requirements; and a two- 
dimensional fit to obtain ipD^X yields (see Fig. 2 
of Ref. jl]). The cross-section thus obtained is model- 
dependent, relying on simulated cc — ► D^X fragmen- 
tation by PYTHIA to establish efficiencies. There has 
been no published update of this remarkable result. 

(Model-independent but since-unpublished Belle re- 
sults were presented to the Quarkonium Working 
Group in 2002, relying instead on counting of D°, D + , 
and all ground-state charmed hadrons under minimal 
cuts. Resolution was similar, with a lower limit on the 
J/ipcc fraction of 0.48 at 95% confidence.) 

A complementary problem is the nature of e + e _ — > 
J/ipX production by processes other than J/t[>cc. 
Such production does seem to occur: there is a con- 
tinuous component in inclusive recoil mass spectra 
(Fig. [5| in excess of background, below the open- 
charm threshold. The mystery is why this component 
- due to the e + e~ — ► if> g process? — should obey 
the cc threshold. A problem for experimental study 
here is the lack of theoretical guidance: predictions 
for the various e + e~ — ► i/j X processes that so pre- 
occupied the original Belle [1\ and BaBar [2] analyses 
have been discredited, but new predictive studies have 
not taken their place. Such work is overdue. 

6. Summary 

The period of fundamental doubt about double 
charmonium production is now over. Questions 
concerning the experimental method have been ad- 
dressed, and the potential for confusion by e + e~ — > 
7*7* — > J/ipX or other more exotic processes has 
been excluded [8|. And the fear that some untrace- 
able mistake had been made is effectively dispelled by 
BaBar 's confirmation [5] of the Belle results |3]. 
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